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On ■ Abstract. We discuss a new completely integrable case of the time-dependent Schrodinger equa- 

I tion in R" with variable coefficients for a modified oscillator, which is dual with respect to the 

' time inversion to a model of the quantum oscillator recently considered by Meiler, Cordero-Soto, 

^ , and Suslov. A second pair of dual Hamiltonians is found in the momentum representation. Our 

d • examples show that in mathematical physics and quantum mechanics a change in the direction of 

I time may require a total change of the system dynamics in order to return the system back to its 

' ' original quantum state. Particular solutions of the corresponding nonlinear Schrodinger equations 

I are obtained. A Hamiltonian structure of the classical integrable problem and its quantization are 

^ also discussed. 

^ : 

c3 , 1. Introduction 

S: 

The Cauchy initial value problem for the Schrodinger equation 

>; t^ = H{t)^ (1.1) 

^ ■ for a certain modified oscillator is explicitly solved in Ref. [ID] for the case of n dimensions in 

ff^ ■ When n = 1 the Hamiltonian considered by Meiler, Cordero-Soto, and Suslov has the form 

O ■ H{t) = - (aat + a^a) + -e^^'a" + -e'^^* {a))\ (1.2) 

00 ' 2 2 2 

O where the creation and annihilation operators are defined as in |28j : 



J3 
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f d \ I f d 
X \ . a = — = - — 



The corresponding time evolution operator is found in [IQ] as an integral operator 

POO 

^ (x, t) = U (t) (x, 0) = / G (x, y, t) ^ (y, 0) dy (1.4) 



with the kernel (Green's function or propagator) given in terms of trigonometric and hyperbolic 
functions as follows 

G{x,y,t) = ^ = (1.5) 

y 27ri (cos t sinh t + sin t cosh t) 



X exp 



2 



y"^) sin t sinh t + 2xy — (x^ + y'^) cos t cosh t 
2i (cos t sinh t + sin t cosh t) 
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It is worth noting that the time evolution operator is known exphcitly only in a few special 
cases. An important example of this source is the forced harmonic oscillator originally considered 
by Richard Feynman in his path integrals approach to the nonrelativistic quantum mechanics [23], 
pi] , [25] , [26] , and [27] ; see also [37] . Since then this problem and its special and limiting cases were 
discussed by many authors; see Refs. [5], [29], [32], [39], [12], [65] for the simple harmonic oscillator 
and Refs. [2], [8], [31], [15], [53] for the particle in a constant external field and references therein. 
Furthermore, in Ref. [20] the time evolution operator for the one-dimensional Schrodinger equation 
(11. ip has been constructed in a general case when the Hamiltonian is an arbitrary quadratic form 
of the operator of coordinate and the operator of linear momentum. In this approach, the above 
mentioned exactly solvable models, including the modified oscillator of [10] , are classified in terms of 
elementary solutions of a certain characteristic equation related to the Riccati differential equation. 

In the present paper we find the time evolution operator for a "dual" time-dependent Schrodinger 
equation of the form 

z^=H{r)ij, r = lsinh(2t) (1.6) 

OT 2 

with another "exotic" Hamiltonian of a modified oscillator given by 

H{t) = ^ (aat + a^a) + ^e~i^'^^^M2r)^2 ^ lgiarctan(2r) (^^t^2 ^ ^^^^^ 

We show that the corresponding propagator can be obtain from expression (II. 5p by interchanging 
the coordinates x ^ y. This implies that these two models are related to each other with respect 
to the inversion of time, which is the main result of this article. 

The paper is organized as follows. In section 2 we derive the propagators for the Hamiltonians 
(ll.2p and (11.70 following the method of [20] — expression (11.50 was obtain in [40j by a totally 
different approach using SU (1, l)-symmetry of the ra- dimensional oscillator wave functions and 
the Meixner-Pollaczek polynomials. Another pair of completely integrable dual Hamiltonians is 
also discussed here. The "hidden" symmetry of quadratic propagators is revealed in section 3. 
The next section is concerned with the complex form of the propagators, which unifies Green's 
functions for several classical models by geometric means. In section 5 we consider the inverses of 
the corresponding time evolution operators and its relation with the inversion of time. A transition 
to the momentum representation in section 6 gives the reader a new insight on the symmetries of the 
quadratic Hamiltonians under consideration together with a set of identities for the corresponding 
time evolution operators. The n-dimensional case is discussed in sections 7 and 8. Particular 
solutions of the corresponding nonlinear Schrodinger equations are constructed in section 9. The 
last section is concerned with the ill-posedness of the Schrodinger equations. Three Appendixes at 
the end of the paper deal with required solutions of a certain type of characteristic equations, a 
quantum Hamiltonian transformation, and a Hamiltonian structure of the characteristic equations 
under consideration, respectively. 

As in [20], [ID] and [61], we are dealing here with solutions of the time-dependent Schrodinger 
equation with variable coefficients. The case of a corresponding diffusion-type equation is investi- 
gated in [02]- These exactly solvable models are of interest in a general treatment of the nonlinear 
evolution equations; see ^, [H], [IB], [H], [22], [SI], [H], [M] and [3], [lU], [II], [12], [13], [IS], 
fl6\ . [33], [44j, [50], [51], [55], [57], [HO] and references therein. They facilitate, for instance, a 
detailed study of problems related to global existence and uniqueness of solutions for the nonlin- 
ear Schrodinger equations with general quadratic Hamiltonians. Moreover, these explicit solutions 
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can also be useful when testing numerical methods of solving the Schrodinger and diffusion-type 
equations with variable coefficients. 

2. Derivation of The Propagators 

The fundamental solution of the time-dependent Schrodinger equation with the quadratic Hamil- 
tonian of the form 



i-^ = -a{t)-^ + b (t) x'lP -tic{t) x-^ + d{t)ij 



in two interesting special cases, namely, 

a = cos^ t, b = sin^ t, c = 2d = sin (2t) (2.2) 

and 

a = cosh^ t, b = sinh^ t, c = 2d = — sinh (2t) , (2-3) 

corresponding to the Hamiltonians (11.21) and (11.71) . respectively, (we give details of the proof in the 
Appendix B), can be found by the method proposed in [20] in the form 



^ = G{x,y,t) = ^. ^^{ait).^+m-y+,it)y^)^ 

A/27ri/i (t) 

where 

^ _ d(ti 

^ ' 4a (t) fi (t) 2a (t) ' 
^ ^ ^ (t) ' dt /i (tf 

/i(t)/i'(t) 2a (0) Jo (/i'(r))^ 
and the function /i (t) satisfies the characteristic equation 

fi" - T {t) fi' + Aa{t)fi = 

with 

T{t) = 2c + 4c/, a (t) = a6 - erf + rf^ + 

a 2 \ a d 



2.1) 



2.4) 

2.5) 
2.6) 
2.7) 

2.8) 
2.9) 



subject to the initial data 

^(0) = 0, /i' (0) = 2a (0) ^ 0. (2.10) 
Equation (12. 5p (more details can be found in [20]) allows us to integrate the familiar Riccati nonlinear 
differential equation emerging when one substitutes (12. 4p into (12. ip . See, for example, [30], ^3] . 
[51], [52], [68] and references therein. A Hamiltonian structure of these characteristic equations is 
discussed in Appendix C. 

In the case (12. 2p . the characteristic equation has a special form of Ince's equation [38] 

/i" + 2tant /i'-2/i = 0. (2.11) 
Two linearly independent solutions are found in |20j : 

Hi = cost cosh t + sin t sinh t = (cost, sinh t) , (2.12) 
fj,2 = cost sinh t + sin t cosh t = (cost, cosh t) (2.13) 
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with the Wronskian W {fJ-i, (^2) — 2cos^t = 2a. Another method of integration of all characteristic 
equations from this section is discussed in the Appendix A; see Table 1 at the end of the paper for 
the sets of fundamental solutions. The second case fl2.3p gives 

/i" -2tanht /i' + 2/i = (2.14) 

and the two linearly independent solutions are [62] 

fi2 = COS t sinh t + sin t cosh t = 14^ (cost, cosh t) , (2.15) 
/i3 = sin t sinh t — cost cosh t = (sin t, cosh t) (2.16) 

with W {fi2^ f^s) = 2cosh^t = 2a. Equation fl2.14p can be obtain from (12. lip as a result of the 
substitution t it. Also, W {fii, fi^) = sin (2t) + sinh (2t) . The common solution of the both 
characteristic equations, namely, 

/i (t) =1^2 = costsinht + sintcosht, (2-17) 

satisfies the initial conditions fl2.10p . 

From fi2.4p - fi2.7p . as a result of elementary calculations, one arrives at the Green function (11.50 
in the case (12.20 and has to interchange there x <-> ?/ in the second case (12. 3p . The reader can see 
some calculation details in section 9, where more general solutions are found in a similar way. The 
next section explains this unusual symmetry between two propagators from a more general point 
of view. 

Two more completely integrable cases of the dual quadratic Hamiltonians occur when 

a = sm^t, b = cos'^t, c = 2d = -sm{2t) (2.18) 

and 

a = sinh^ t, b = cosh^ t, c = 2d = sinh (2t) . (2.19) 
The corresponding characteristic equations are 

/i"-2cott/i'-2/i = (2.20) 

and 

/i"-2cotht/i' + 2/i = 0, (2.21) 
respectively, with a common solution 

yU (t) = = sin t cosh t — cos t sinh t = W (sin t, sinh t) (2.22) 

such that /i (0) = /i' (0) = fi" (0) = and (0) = 4. Once again, from fl2.4p - fl2.7p one arrives at 
the Green function 

G{x,y,t) = ^ = (2.23) 
^y 2ni (sin t cosh t — cos t sinh t) 

[x^ + u^) cos t cosh t — 2xy + (x^ — u^) sin t sinh t 

X exp ' 



2i (cos t sinh t — sin t cosh t) 

in the case f l2.18p and has to interchange there x ^ y in the second case (12.190 . The corresponding 
asymptotic formula takes the form 
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as e = 0+. We will show in section 6, see Eqs. (16. ip and (16.81) . that our cases (I2.2p - (I2.3I) and 

fl2.18p - (l2.19p are related to each other by means of the Fourier transform. 

We have considered some elementary solutions of the characteristic equation (12. Sp . which are of 
interest in this paper. Generalizations to the forced modified oscillators are obvious; see Ref. |40j . 
More complicated cases may include special functions, like Bessel, hypergeometric or elliptic func- 
tions tlj, n [32], [63], and [68]. 



3. On a "Hidden" Symmetry of Quadratic Propagators 

Here we shall elaborate on the symmetry of propagators with respect to the substitution x ^ y. 
Lemma 1. Consider two time- dependent Schrodinger equations with quadratic Hamiltonians 

i^ = -akit)^ + hit)x''ij-i(^Ckit)x^ + dk{t)ij'^ {k = l,2), (3.1) 

where ci — 2di = C2 — 2^2 = £ (t) and dk (0) = 0. Suppose that the initial value problems for 
corresponding characteristic equations 

fi" -Tk{t)fi' + Aak{t)fi = 0, /i(0)=0, /x' (0) = 2afc (0) ^ (3.2) 

with 

Tk (t) = — - 2cfc + 44, CTk (t) = akbk - Ckdk + dl + ^ (^-^] (3.3) 
flfc 2 \ak dkj 

have a joint solution fi (t) and, in addition, the following relations hold 

4 {a^b^ - c,d, + df) = ^"^"^^^ ~ " = ^ («2&2 - C24 + dj) , (3.4) 

where h (t) = exp ^— J e (r) dr^ . Then the corresponding fundamental solutions 

= Gk {x, y, t) = ^^{Mt>'+Mty-y+l,{t)y') (3,5) 

^y2Tcifj, (t) 

possess the following symmetry 

a (t) = «i (t) = 72 (t) , 7 (t) = 7i (t) = «2 (t) , /3{t)= (3, (t) = [3, (t) (3.6) 

and 

Gi{x,y,t) = G2iy,x,t) . (3.7) 
This property holds for a single Schrodinger equation under the single hypothesis Ili3.4\ )- 

Indeed, according to Ref. [2U] . 

f3{t) = (3,{t)=P2it) = -^ (3.8) 

in the case of a joint solution fi (t) of two characteristic equations. In view of the structure of 
propagators for general quadratic Hamiltonians found in [20], the symmetry under consideration 
holds if we have 

1 a' di da h'^ ^ ,^ 

4ai jj la\ dt 
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and 

1 = -. 7^, 3r + «i— = 0, 3.10 

4a2 /i 2a2 at jj,"^ 

simultaneously. Excluding a from (13. 9p . one gets 

/ - ^/z' + 2rf, f ^ - ^ V = (3.11) 

Comparison with the characteristic equation results in the first condition in (13.41) . The case of 7, 
which gives the second condition, is similar. This completes the proof. 

A few examples are in order. When a = 1/2, 6 = c = (i = 0, and /i" = 0, = t, one gets 

as the free particle propagator [27] with an obvious symmetry under consideration. Our criteria 
(13. 4p . namely, Aa? = (/x')^ , stands. 

The simple harmonic oscillator with a = b = 1/2, c = d = and fi" + fi = 0, fi = sint has the 
familiar propagator of the form 

Gix,y,t) = ^=L== exp(-^ {{x^ + y^)cost-2xy)] , (3.13) 
V27rzsmt V2smt V 

which is studied in detail at [5], [29], [32], [39], [12], [65]. (For an extension to the case of the forced 
harmonic oscillator including an extra velocity-dependent term and a time- dependent frequency, 
see [23], [21], [2Z] and [37].) Our condition (13.41) takes the form of the trigonometric identity 

4„,.l;^, (3.14) 

which confirms the symmetry of the propagator. 

For the quantum damped oscillator [21] a = b = ujq/2, c = d = —A and 



G (x, y,t) = \ exp [[x + y j cos ut — 2xyj 

^' 2muJosmLjt \2ci;osma;t 

xexp(|^(x^-y^)) (3.15) 
with uj = a/cUq — > and /i = (cuq/uj) e~^^ sin tot. The criterion 

4„6=iMl^-2.^^:, (3.16) 

where e = c — 2d = X and h = e~'^*, holds. But here ci (0) = — A 7^ and a more detailed analysis of 
asymptotics gives an extra antisymmetric term in the propagator above; see [21] for more details. 

The case of Hamiltonians (11.21) and (11.71) corresponds to 

ai = cos^t, bi = sin^t, ci = 2di = sin {2t) (3-17) 

and 

a2 = cosh^ t, 62 = sinh^ t, C2 = 2^2 = — sinh (2t) , (3.18) 
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respectively. Our criteria (13 ■4p are satisfied in view of an obvious identity 

4aia2 = {fi'f . (3.19) 

Tlie cliaracteristic function is given by (12.171) . Tliis explains the propagator symmetry found in the 
previous section. 

Our last dual pair of quadratic Hamiltonians has the following coefficients 

ai = sin^t, 6i = cos^t, Ci = 2(ii = — sin (2t) (3.20) 

and 

a2 = sinh^t, 62 = cosh^t, C2 = 2^2 = sinh (2t) . (3-21) 

The criteria (13. 4p are satisfied in view of the identity (I3.19P with the characteristic function (12.221) 
and, therefore, the propagator (I2.23P obeys the symmetry under the substitution x y. 

Remark 1. A simple relation 

/i T1-T2 

which is valid for a joint solution of two characteristic equations, can he used in our criteria l[3.4\}- 



Although we have formulated the hypotheses of our lemma for the Green functions only, it can be 
applied to solutions with regular initial data. For instance, a pair of characteristic equations (12.111) 
and (12.211) has a joint solution given by (I2.12p . which does not satisfy initial conditions required for 
the Green functions. The coefficients of the corresponding quadratic Hamiltonians are 

Oi = cos^t, hi = sin^t, Ci = 2di = sin (2t) (3.23) 

and 

a2 = sinh^ t, 62 = cosh^ t, C2 = 2^2 = sinh (2t) . (3.24) 
The criteria (13. 4p are satisfied in view of the identity (I3.19P and the particular solution 

iIj = K{x, y, t) = ^ = (3.25) 

2tt (cos t cosh t + sin t sinh t) 



{x^ + y^) sin t cosh t — 2xy — {x^ — y^) cos t sinh t 

X exp ' 



2i (cos t cosh t + sin t sinh t) 

obeys the symmetry under the substitution x ^ y. The initial condition is the standing wave 
K{x,y,0) = e^^Vv^. 

In a similar fashion, the characteristic equations (I2.20p and (I2.14p have a common solution 

fi = — /ig = cost cosh t — sin t sinh t. (3.26) 
The coefficients of the corresponding Hamiltonians are 

ai = sin^ t, hi = cos^ t, ci = 2di = — sin (2t) (3.27) 

and 

02 = cosh^ t, 62 = sinh^ t, C2 = 2^2 = - sinh (2t) . (3.28) 
The criteria (13. 4p are satisfied once again and the particular solution is given by 

i(j = K {x,y,t) = ^ (3.29) 

y 2n (cos t cosh t — sin t sinh t) 
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[x + y ) sin t cosh t + 2xy + (x — y ) cos t sinh t 
2i (cos t cosh t — sin t sinh t) 



with i^' (x, 0) = e ^'^^Z -\/27r. We shall discuss in section 6 how these solutions are related to the 
corresponding time evolution operators. 

4. Complex Form of The Propagators 

It is worth noting that the propagator in (11.51) can be rewritten in terms of the Wronskians of 
trigonometric and hypergeometric functions as 

\J iTxiW (cost, cosh rj 



X exp 



W (sin t, cosh t) + 2xy — W (cos t, sinh t) 



2iW (cost, cosht) 
This simply means that our propagator has the following structure 



G = exp (^^^' + ^^^-^^^^^-^^^') , (4.2) 

where the coefficients are solutions of the system of linear equations 

ci cos t + C2 cosh t = sint, (4.3) 
—Ci sin t + C2 sinh t = cost, 

C3 cos t + C4 cosh t = sinh t, 
—C3 sin t + C4 sinh t = cosht 

obtained by Cramer's rule. A complex form of this system is 

ciz* + C2C = iz* , C32;* + C4C = <*, (4.4) 

where we introduce two complex variables 

^ = cost + i sin t, = cosh t + i sinh t (4.5) 

and use the star for complex conjugate. Taking the complex conjugate of the system (14. 4p . which 
has the real-valued solutions, and using Cramer's rule once again, one gets 

C3 = —r-^ -qrT' C4 



i{zC-z*Cy i{zC-z*C)' 

As a result, we obtain a compact symmetric expression of the propagator (11. 5p as a function of two 
complex variables 

G{x,y,t) = G{x,y,zX)= , , ] (4.7) 

zC + -2*C*) - 4xy + {zC + -2*0 



X exp 

This function takes a familiar form 



2 (^*C - zQ 
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z = cost+ismt 





Figure 1. Two synchronized "clocks", or contours in complex Euclidean z and 
pseudo-Euclidean ( "time" planes, corresponding to modified oscillators. 



G 



1 



: exp 



(x^ — y^) X2X4 + 2xy — {x^ + y"^) Xix-^ 
2i (xiXi + X2X3) 



(4. 



in a real- valued four- dimensional vector space, if we set z = Xi + ix2 and ( = x^ + ix^ with x[ = —X2, 



2 — I5 3 



X4, x'^ = X3, and solve the following initial value problem 



x'l + Xi 


= 0, 


xi (0) 


= 1, 


x\ (0) 


= 0, 


(4.9) 


X2 + X2 


= 0, 


X2 (0) 


= 0, 


x'2 (0) 


= 1, 




X3 — X3 


= 0, 


X3 (0) 


= 1, 


4 (0) 


= 0, 




Ob ^ X4 


= 0, 


Xa (0) 


= 0, 


< (0) 


= 1, 





the solution of which can be interpreted as a trajectory of a classical "particle" moving in this space; 

cf. (USD. 

It is worth noting that our propagators expression (14.71) . extended to a function of two indepen- 
dent complex variables z and C, allows us to unify several exactly solvable quantum mechanical 
models in geometrical terms, namely, by choosing different contours, with certain "synchronized" 
parametrization, in the pair of complex "time" planes under consideration. Indeed, the free particle 
propagator (I3.12p appears in this way when one chooses z = 1 and = 1 + it. The simple harmonic 
oscillator propagator (I3.13P corresponds to the case z = 1 and ( = e**. As we have seen in this 
section, the propagator (II. 5p is also a special case of (14.71) . This is why we may refer to the Hamil- 
tonians under consideration as the ones of modified oscillators. By a vague analogy with the special 
theory of relativity, one may also say that in this case there are two synchronized "clocks" , namely, 
the two contour parameterized by (14. 5p . one in Euclidean and another one in the pseudo- Euclidean 
two dimensional spaces, respectively, which geometrically describes a time evolution for the Hamil- 
tonians of modified oscillators; see Figure 1. This idea of introducing a geometric structure of time 
in the problem under consideration may be useful for other types of evolutionary equations. 



In a similar fashion, our new propagator (12.230 can be rewritten in terms of the Wronskians as 

1 



G{x,y,t) 



^j2TxiW (sint, sinht) 



(4.10) 



10 



RICARDO CORDERO-SOTO AND SERGEI K. SUSLOV 




u-Kt) 



Figure 2. The time evolution operator and its inverse. 



X exp 

The corresponding complex form is 



W (cos t, sinh t) — 2xy — W (sin t, cosh t) 



G{x,y,t) = G{x,y,z,C) 



2iW (sinh t, sin t) 
1 



X exp 

We leave the details to the reader. 



(zC + zX) x^-4xy + {z( + z*C)y^ 
2{z%-zC) 



(4 



5. The Inverse Operator and Time Reversal 

We follow the method suggested in ^61j for general quadratic Hamiltonians with somewhat 
ferent details. The left inverse of the integral operator defined by fll.4p - fll.5l) . namely, 



u (t) {x) 



G{x,y,t) i){y) dy, 



is 



U-'{t)xix) = / G{y,x,-t) xiy) dy 

J — oo 

in view of U^^ = f/^; see Figure 2. Indeed, when s < t, by the Fubini theorem 

POO 

U-\s){U{t)ij) = U'\s)x= G{z,x,-s) x{z) dz 



G{z,x,-s) yj G{z,y,t) i){y) dy ] dz 
G{x,y,s,t) ^piy) dy. 



Here 



G{x,y,s,t) 



G {z,x, —s) G {z,y,t) dz 



TIME INVERSION FOR OSCILLATORS 11 

^ / ^i{{a(t)-a(s))z^ + (my-m^>) 

27r^/i(s)/i(t) 

1 



y^Anifi (s) fi (t) {a (s) — a (t)) 

\f3{t)y-(3 (s) xf -A {a it) -a (.)) (7 it) - 7 {s) x') 



X exp 



Ai {a (t) - a (s)) 

by the familiar Gaussian integral [6], [58] and [56]: 

00 



In view of (12.61) . 



/i (s) jJ'it) J fi (s) fi (t) {a (t) — a (s)) 



{f3{t)y-f3{s)xf = (^--^]'= ' , (5.6) 



and a singular part of (15. 4p becomes 



1 {f3{t)y-f3{s)x) 
exp 



2 



v/47r^/i (s) /X (t) (a(s)-a(t)) I 4^ (a (t) - a (s)) 



1 

exp 



Amiyi, (s) [i it) {a (s) — a it)) 



^4i/i (s) (t) (a (t) — « (s)) ^ 



Thus, in the limit s ^ t , one can obtain formally the identity operator in the right hand side of 
(15. 3p . The leave the details to the reader. 

On the other hand, the integral operator in (ll.4l) - (ll.5p . namely, 

X (x) = ^=i== r e<"(*)^^+^W^^+-W^^) ^ iy) dy (5.7) 
l-ni^i [tj J-oo 

is essentially the Fourier transform and its inverse is given by 

^ {y) = ^=i= r e-<"W^^+^W^^-^-W^^) X {x) dy (5.8) 

in correspondence with our definition (15. 2p in view of (12. 6p . 

The Schrodinger equation (11.10 retains the same form if we replace t in it by —t and, at the 
same time, take complex conjugate provided that {H {—t) (f)* = H(t)(f*. The last relation holds 
for both Hamiltonians (ll.2p and (ll.7p . Hence the function x {x, t) = ip* (x, — t) does satisfy the same 
equation as the original wave function ip (x, t) . This property is usually known as the symmetry 
with respect to time inversion (time reversal) in quantum mechanics [29j, [36j, |l2], [69]. This fact 
is obvious from a general solution given by (11.40 - 01.50 for our Hamiltonians. 
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On the other hand, by the definition, the (left) inverse U ^ {t) of the time evolution operator 
U {t) returns the system to its initial quantum state: 

^{x,t) = U {t)ij{x,0) , (5.9) 
U-^ (t) ifj {x, t) = U-^ (t) {U (t) ip {x, 0)) = ifj (x, 0) . (5.10) 

Our analysis of two oscillator models under consideration shows that this may be related to the 
reversal of time in the following manner. The left inverse of the time evolution operator (11. 4p for 
the Schrodinger equation (11.11) with the original Hamiltonian of a modified oscillator (11.21) can be 
obtained by the time inversion t — — t in the evolution operator corresponding to the new "dual" 
Hamiltonian (II. 7p (and vise versa). The same is true for the second pair of dual Hamiltonians. 
More details will be given in section 6. This is an example of a situation in mathematical physics 
and quantum mechanics when a change in the direction of time may require a total change of 
the system dynamics in order to return the system back to its original quantum state. Moreover, 
moving backward in time the system will repeat the same quantum states only when 

^{x,t-s) = U{t-s)'ip{x,0) = U-^{s)U{t)ij{x,0), 0<s<t, (5.11) 

which is equivalent to the semi-group property 

U{s)U{t-s) = U{t) (5.12) 

for the time evolution operator. This seems not true for propagators (II. 5p and (I2.23p . 



6. The Momentum Representation 



The time-dependent Schrodinger equation (12. ip can be rewritten in terms of the operator of 
coordinate x and the operator of linear momentum = i~^d/dx as follows 

oih 

i-^ = {a{t)pl + b{t)x^ + d{t){xp.^+p^x))^ (6.1) 

with c = 2d. The corresponding quadratic Hamiltonian 

H = a{t)pl + b (t) x'^ + d (t) {xp^ + p^x) (6.2) 

obeys a special symmetry, namely, it formally preserves this structure under the permutation x ^ 
Px- This fact is well-known for the simple harmonic oscillator [29], [36], [12] . 

In order to interchange the coordinate and momentum operators in quantum mechanics one 
switches between the coordinate and momentum representations by means of the Fourier transform 

1 r°° ■ 

^{x) = -= e'-yx {y) dy = F [x] (6.3) 

V ZTT J -oo 

and its inverse 

X (y) = -= / e-^y^J (x) dx = [i:] . (6.4) 

Indeed, the familiar properties 

Px-ip = PxF [x] = F [yx] , xij = xF [x] = -F [pyx] (6.5) 

imply 

pli; = F [y\] , x'ij = F [pjx] (6-6) 
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and 

{xp^ + p^x) ij = -F [{pyy + ypy) x] ■ (6.7) 

Therefore, 

Hijj = [apl + bx'^ + d {xpx + Pxx)) F [x] 
= F [ [hpl + ay^ -d {ypy + Pyy)) x] 
by the hnearity of the Fourier transform. In view of 



dt 



the Schrodinger equation (16. ip takes the form 



= {b (t) Pl + a (t) y^-d (t) {ypy + pyy)) x (6.8) 
with a b and d —>■ — d in the momentum representation. 

This property finally reveals that our quadratic Hamiltonians ( \2.2\\ and (12.181) . similarly (12.31) 
and (12.191) . corresponds to the same Schrodinger equation written in the coordinate and momentum 
representations, respectively. Thus, in section 2, we have solved the Cauchy initial value problem 
for modified oscillators both in the coordinate and momentum representations. 

In this paper the creation and annihilation operators are defined by 

^ _ px + ix _ f d 



with the familiar commutator [a, a^"] = aa^ — a^a = 1 [28j. One can see that 

a^^p = a^F [x] = F [ittyx] , (6.11) 
O'l'^ = 4^ [x] = F [-ialx] , 

or 

Qx — > idy, a\. (6.12) 
under the Fourier transform. This observation will be important in the next section. 

Finally one can summarize all results on solution of the Cauchy initial value problems for the 
modified oscillator under consideration in a form of the commutative evolution diagram on Figure 3. 
We denote 



U{t)^{x) = / Gu{x,y,t) ^{y) dy, (6.13) 

J — oo 

/oo 

Kuix,y,t) ij{y) dy, (6.14) 
'OO 

/oo 

Gv{x,y,t) ij{y) dy, (6.15) 
-oo 

POO 

L{t)^{x) = / Kvix,y,t) ^Piv) dy. (6.16) 
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Figure 3. The commutative evolution diagram in R. 



The kernels of these integral operators are defined as follows. Here Gu {x, y, t) and Gy (x, y, t) are 
the Green functions in f ll.Sp and (12.231) . respectively. The kernels Ku{x,y,t) and Kv{x,y,t) are 
given by (13.251) and (I3.29p . respectively. The following operator identities hold 



U{t) = K it) F-^ = FL (t) = FV (t) F-\ (6.17) 

V{t) = L (t) F = F^^K it) = F-^U it) F, (6.18) 

U-^ (t) = FK-^ (t) = L-^ (t) F~^ = FV-^ (t) F^\ (6.19) 

y-i (t) = F-^L-^ (t) = K-^ (t) F = F~^U-^ (t) F, (6.20) 

K{t)= FL{t)F, L{t) = F~^K{t)F-\ (6.21) 

R-^ (t) = F-^L-^ (t) F-\ L-^ (t) = FR-^ (t) F. (6.22) 



Here F and F ^ are the operators of Fourier transform and its inverse, respectively, which relate 
the wave functions in the coordinate and momentum representations 

i^ = F[x], X = F-'[^] 

at any given moment of time thus representing the vertical arrows at our diagram. The time 
evolution operators U (t) , V (t) and their inverses (t) , {t) correspond to the horizontal 
arrows in the coordinate and momentum representations, respectively. They obey the symmetry 
with respect to the time reversal, which has been discussed in section 5. 

In order to discuss the diagonal arrows of the time evolution diagram on Figure 3, we have to 
go back, say, to the particular solution (I3.25p . A more general solution of the Schrodinger equation 
(I2.ip - (I2.2I) can be obtained by the superposition principle in the form 

/oo 
Kuix,y,t) xiy,0) dy, (6.23) 
-oo 

where x is an arbitrary function, independent of time, such that the integral converges and one can 
interchange the differentiation and integration. In view of the continuity of the kernel at t = 0, we 

get 

^ (x, 0) = e'^y X (y, 0) dy, (6.24) 
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which simply relates the initial data in the coordinate and momentum representations. Then 
solution of the initial value problem is given by the inverse of the Fourier transform 



X (y, 0) = -= / e-'^y ij (x, 0) dx (6.25) 
V Zir J-oo 

followed by the back substitution of this expression into fl6.23p . This implies the above factorization 
U (t) = K (t) of the corresponding time evolution operator; see (16.171) . The Green function 
(II. 5p can be derived as 

1 

Gu {x, y, t) = -j= / Ku {x, z, t) e"*^" dz (6.26) 

with the help of the integral (15. 5p . The second equation, U {t) = FL{t) , is related to following 
integral 

1 



Gu (x, y,t) = ^ I e"^ Ky {z, y, t) dz. (6.27) 
V 27r 



CXD 



The meaning of the operator L (t) , represented by another diagonal arrow on the time evolution 
diagram, is established in a similar fashion. One can see that the relation V [t) = L (t) F in (16.181) 
follows from the elementary integral 

1 r°° 

Gv (x, y, t) = -j= / Kv (x, z, t) e'^' dz (6.28) 
and K (t) = FV (t) corresponds to 

1 ■ 

Ku (x, y,t) = ^ e^"^ Gy {z, y, t) dz. (6.29) 

This proves (I6.17p - (l6.18p . The inverses K^^ {t) and L^^ (t) are found, for instance, by the inverse 
of Fourier transform similar to (15.70 - 05.80 . They are not directly related to the reversal of time. 
We leave further details about the structure of the commutative diagram on Figure 3 to the reader. 

7. The Case of n- Dimensions 

In the case of i?" with an arbitrary number of dimensions, the Schrodinger equation for a modified 
oscillator (11.10 with the original Hamiltonian 



s=l s=l s=l 

considered by Meiler, Cordero-Soto, and Suslov [IQ], has the Green function of the form 

n 

Gt{x,x') = l[Gt{xs,x'^) (7.2) 

s=l 



27ri (cos t sinh t + sin t cosh t) 

— x'"^) sin t sinh t + 2x ■ x' — ix"^ + x'"^) cos t cosh t 



X exp 



2i (cos t sinh t + sin t cosh t) 
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Solution of the Cauchy initial value problem can be written as 

ilj{x,t)= Gtix,x') iIj{x',0) dx', (7.3) 

where dv' = dx' = dx'i ■ ... ■ dx'^. The propagator expansion in the hyperspherical harmonics is given 
by 

Gt {x, x') = Yku m Yk. m Qf (r, r') (7.4) 



with 



Kv 



g-i7r(2ii'+n)/4 [rr'^^ 



^''^"'"'^ 2^+-/2-ir (ir + n/2) (costsinht + sin t cosh t)^+"/^ ^^'^^ 



X exp 



(r^ _l_ (^^'~) ^ cost cosh t — (r^ — (r') ) sintsinht 



2 (cos t sinh t + sin t cosh t) 



JM = ^^T^.oF,( ; --] (7.7) 



K + n/2 ' 4 (cost sinh t + sin t cosh t)^ 

Here Ykv {^) are the hyperspherical harmonics constructed by the given tree T in the graphical 
approach of Vilenkin, Kuznetsov and Smorodinskii ^6], the integer K corresponds to the constant 
of separation of the variables at the root of T (denoted by K due to the tradition of the method 
of i^'-harmonics in nuclear physics [22]) and v = {/i,/2,--- Jp} is the set of all other subscripts 
corresponding to the remaining vertexes of the binary tree T. These formulas imply the familiar 
expansion of a plane wave in i?" in terms of the hyperspherical harmonics 

— E (^) Yk. m {rr') , (7.6) 

where 

r(/i + i) °^ + 1 ' 4 

is the Bessel function. See [ID] and references therein for more details. It is worth noting that 
the Green function (11.51) was originally found by Meiler, Cordero-Soto, and Suslov as the special 
case n = 1 of the expansion fl7.4l) - fl7.5l) . The dynamical SU (1,1) symmetry of the harmonic 
oscillator wave functions, Bargmann's functions for the discrete positive series of the irreducible 
representations of this group, the Fourier integral of a weighted product of the Meixner-PoUaczek 
polynomials, a Hankel-type integral transform and the hyperspherical harmonics were utilized in 
order to derive the n-dimensional Green function. 

Our results show that the "dual" Schrodinger equation (11.61) with a new Hamiltonian of the form 

^ n n ^ n 

^ (^) = ^ E + + ^e-^---(^-) («^)' + ^e^^--(^^) E (7-8) 

S = l S = l 3 = 1 

has the propagator that is almost identical to (17. 2p but with x ^ x' . Indeed, in the case of n- 
dimensions one has 

n 

H{r) = YHs{r), (7.9) 

s=l 
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where we denote 



If one chooses 



'^s = i>s{xs,t) = Gtix'^Xs) (7.11) 



1 \ 1/2 



27ri (cos t sinh t + sin t cosh t) 

x'/ — x"^) sin t sinh t + 2x'^Xs — {xf + x^) cos t cosh t 



X exp 



with 



and denotes 



then 



and 



As a result, 



2i (cos t sinh t + sin t cosh t) 

dib 1 
2^ = i7,(r)^„ t = -sinh(2r) (7.12) 



^ = ll'^k = llGt {xl Xk) = Gt {x', x) , (7.13) 

k=l k=l 

s=l ^ ' fc^s 

n 

s=l k^s 



d 



and Eq. (11.61) for the n-dimensional propagator is satisfied. For the initial data, formally, 

n n 

lim Gt {x', x) = Y\ 1™ {x'f,, Xk) = W^ Wk " ^k) = S {x — x') , (7.17) 



fc=l k=l 



where 6 {x) is the Dirac delta function in ii". Further details are left to the reader. 

The n-dimensional version of the Hamiltonian corresponding to the coefficients (12.181) is given by 



with the propagator 



2 ^ ' 2 

s=l s=l s=l 



I \ «/2 



Gt(x,x') = — (7.19) 

\ 2m (sm t cosh t — cos t smh t) ) 

( [x^ + a;'^) cos t cosh t — 2x-x'^ [x^ — x'"^) sin t sinh t 

X exp ^ ^ — - — - — 

\ 2% (cos t sinh t — sin t cosh t) 
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The dual counterpart of this Hamiltonian with respect to time reversal has the form 

s=l s=l s=l 

and one has to interchange x x' in (17.191) in order to obtain the corresponding Green function. It 
is worth noting that the Hamiltonians (17.11) and fl7.18p (respectively, (17.81) and (I7.20p ) are transform- 
ing into each other under the substitution — > ias-, a\ — >■ — ^oj, which preserves the commutation 
relations of the creation and annihilation operators. As we have seen in the previous section this 
property is related to solving the problem in the coordinate and momentum representations. 

Combining all four cases together, one may summarize that two Hamiltonians, 

s=l s=l s=l 

and their duals with respect to the time reversal. 



s=l s=l s=l 

with T = I sinh (2t) , have the following Green functions: 

\ 2m (sm t cosh t ± cos t smh r j / 

/ ± (x"^ — x'"^) sin t sinh t + 2a:; ■ as' — (a;^ + £c'^) cos t cosh t 

X exp — ^ ; ^ — - — r — 

\ 2i (sin t cosh t ± cos t sinh t) 

This expression is valid for two Hamiltonians (17.211) . respectively. One has to interchange x ^ x' 
for the case of the dual Hamiltonians (17.221) . 

In a similar fashion, the n-dimensional form of the kernels (13.251) and (13.291) is 

Kf {x, x') = — (7.24) 

\ 27r (cos t cosh t ± sin t sinh t) J 

{x'^ + a?'^) sin t cosh t ^ 2x ■ x' ^ (a?^ — x''^) cos t sinh t 

X exp ' 



2i (cos t cosh t ± sin t sinh t) 



and 



GH^.^') = TT^TI / i^,^(^,a;") e^^-'-"rfx" (7.25) 
(27r) ' Jh" 

1 f g±ia=.a=" ^a." a;') dx". 



{2ny Jr- 



Denoting 



U±{t)i){x) = I Gf{x,x') i){x') dx', (7.26) 
K±{t)i;{x) = [ Kf{x,x') i){x') dx', 
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Ulit) 



Figure 4. The commutative evolution diagram in i?". 



and 



1 



(27r)"/' JR- 

one arrives at the commutative evolution diagram in IV on Figure 4. The corresponding relations 
are 

f/± it) = it) = F^K^ (t) = F^U^ (t) F^, (7.28) 

U^' (t) = F^Kg' (t) = (t) F^ = F^U-' (t) F^. 

We leave further details to the reader. 

A certain time-dependent Schrodinger equation with variable coefficients was considered in [10] 
in a pure algebraic manner in connection with representations of the group SU (1, 1) in an abstract 
Hilbert space. Our Hamiltonians (17.211) and (I7.22p belong to the same class thus providing new 
explicit realizations of this model in addition to several cases already discussed by Meiler, Cordero- 
Soto, and Suslov. 



if) {x') dx' 



(7.27) 



8. ElGENFUNCTION EXPANSIONS 



The normalized wave functions of the n-dimensional harmonic oscillator 
have the form 



1 

= 9 5Z 



s=l 



(8.1) 



^ {x) = ^NKu (r, Q) = Yk, {Q) Rmk (r) , (8.2) 
where Ykv {^) are the hyperspherical harmonics associated with a binary tree T, the integer number 
K corresponds to the constant of separation of the variables at the root of T and i/ = {/i, ^2, ••• , ^p} 
is the set of all other subscripts corresponding to the remaining vertexes of the binary tree T; see 
[IS] , [59] , [67] for a graphical approach of Vilenkin, Kuznetsov and Smorodinskii to the theory of 
spherical harmonics. The radial functions are given by 



R 



NK 



(r) 



2[{N -K) /2]! 
r [(AT + X + n) /2] 



exp (-rV2) L 



K jK+n/2-l 
(N-K)/2 



(8.3) 
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where are the Laguerre polynomials. The corresponding energy levels are 

E = En = N + n/2, (N - K) /2 = k = 0,1,2, ... (8.4) 

and we can use the SU (1, l)-notation for the wave function as follows 

ijjj^^^y (x) = "^NKu (r, n) = (fi) Rnk (r) , (8.5) 

where the new quantum numbers are given by j = i^/2+n/4 — 1 and m = A^/2+n/4 with m = j + 
j + 2, ... . The inequality m > j + 1 holds because of the quantization rule (18. 4p . which gives N = K, 
K + 2, K + 4, ... . See [10], [IBJ and [59j for more details on the group theoretical properties of the 
n-dimensional harmonic oscillator wave functions. 

The Cauchy initial value problem for the Schrodinger equation (II. ip with the Hamiltonian of a 
modified oscillator (17. ip has also the eigenfunction expansion form of the solution |40j : 

oo 

j{u} m=j+l 

with the time dependent coefficients 

oo „ 

Cm (t) = Yl '"^'"^^'m m / ^ r^m'lu} i^') ^ (^', 0) dv' (8.7) 

m'=j+l -^^^ 

given in terms of the Bargmann functions [1], [16] and [67] 



= r i2j + 2) V (m-j-l)!(m--j-l)! l^^"^ 2 J l'""^ 2 



m+m 



^ ^^^l 2j + 2 ' -sinh^(^/2)J- ^^-^^ 

Choosing the initial data in (17. 3p and (18. 60 - 08. 7p as ip {x,0) = 6 {x — x') , we arrive at the eigen- 
function expansion for the Green function 



oo 



G,{x,x') = J2 E e-'^-U"^'-"^vi„^{2t) ^^^^,^{x)rjm'iu}{^'), (8.9) 

m,m' =j+l 

where by (I7.19P the following symmetry property holds 

Gt{x,x') = G%{x,x'). (8.10) 

In this paper we have found solution of the Cauchy initial value problem for the new Hamiltonian 
(17. 8p in an integral form 

ilj{x,t)= Gtix',x) '^{x',0) dx'. (8.11) 
In view of (l8.9p -( !8rT0l) . the eigenfunction expansion of this solution is given by 

oo 

^M = J2 Yl ^i-M ' (8.12) 

m=j+l 
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where 

oo „ 

(t) = E (-0"' ^"'^"^'^ (^^'™ (-2^))* / „ (^') ^ (^', 0) cit;'. (8.13) 

m'=j+l 

This expansion is in agreement with the unitary infinite matrix of the inverse operator in the basis 
of the harmonic oscillator wave functions; see section 5. 

The cases of the Hamiltonian (17.181) and its dual (17.201) can be investigated by taking the Fourier 
transform of the expansions (I8.6l) - (l8.7p and (I8.12p - (l8.13p . respectively. The corresponding trans- 
formations of the oscillator wave functions are 

t^^'^NK. (x) = / e^^'^-^' ^^Ku ix') dx'. (8.14) 

(27r) ' Jr" 

This can be evaluated in hyperspherical coordinates with the help of expansion (]7.6p -( j7!Tj) . or 
by adding the St/ (1, l)-momenta according to the tree T |16] and using linearity of the Fourier 
transform. One can use 

e--' = (27r)"/' Yl ^Ku (^) Yk. m 5_i (r, r') (8.15) 

Ku 

with 

{tt'^^ ( — (rr')^ \ 

and 

/■oo 
^0 

as a special case of Eqs. (7.3) and (7.6) of Ref. [40j together with the orthogonality property of 
hyperspherical harmonics. We leave further details to the reader. 

9. Particular Solutions of The Nonlinear Schrodinger Equations 

The method of solving the equation (12.11) is extended in |20j to the nonlinear Schrodinger equation 
of the form 



= -a (t) ^ + & it) x'^ - t (t) x^ + d (t) ^j+h (t) \^P\'' s>0. (9.1) 

We elaborate first on two cases (12.21) and (12.31) . A particular solution takes the form 

^ = (x, t) = Kh (x, y, t) = -^L= e<°W-'+/3(*)-J/+7W2/^+«W) ^ = constant, (9.2) 

where equations (I2.5p - (l2.7p hold and, in addition, 

dt^ hit) , , /"* h (r) , 

«:(t) = «:(0)- / ^j\dT, (9.3) 



dt /i^t) ' 7o ij) 

provided that the integral converges. 

In the first case (12.20 . by the superposition principle, the general solution of the characteristic 
equation (12.110 has the form 

/i = Ci/ii {t) + C2/i2 {t) (9-4) 
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COS t (ci cosh t + C2 sinh t) + sin t (ci sinh t + C2 cosh t) 
e* sin t + - + ^ e"* cos t ' 



V 4; ^2 V 4, 

with fi' = 2 cos t (ci sinh t + C2 cosh t) and 

/x(0) = ci, /i'(0) = 2c2. (9.5) 

Then 

, , cos t (ci sinh t + C2 cosh t) — sin t (ci cosh t + C2 sinh t) 

Ot [ty] = (9.0) 

2 (cos t (ci cosh t + C2 sinh t) + sin t (ci sinh t + C2 cosh t)) ' 

/jm = ^i^W ^ (9 7) 

cos t (ci cosh t + C2 sinh t) + sin t (ci sinh t + C2 cosh t) ' 

and 

(0) (cos t sinh t + sin t cosh t) 

T (^) ~ T (^) — (^•^) 

2 (cos t (ci cosh t + C2 sinh t) + sin t (ci sinh t + C2 cosh t) ) 

as a result of elementary but somewhat tedious calculations. The first two equations follow directly 
from (12. 5p and a constant multiple of the first equation (12. 6p . respectively. One should use 

^ + a (t) = 0, (9.9) 



see [20], integration by parts as in (12.71) . and an elementary integral 

dt sinh t 



Ci sinh t + C2 cosh t) C2 (ci sinh t + C2 cosh t) 



+ C (9.10) 



in order to derive (19.81) . 

Two special cases are as follows. The original propagator (11.51) appears in the limit ci ^ 
when /5 (0) = — (ci)~^ and 7(0) = (2ciC2)~^ . The solution with the standing wave initial data 
ip (x, 0) = e"'^ found in [20j corresponds to Ci = 1 and C2 = 0. 

Equation (19.31) can be explicitly integrated in some special cases, say, when h (t) = Xfi' (t) : 

K (0) - ifi^-' (t) - n^-' (0)) , when s ^ 1, 

^it) = { ~) .... (9-11) 

/€(0)-Aln( ^ ) , whens = l. 

V/i(0)y 

Here /i (0) 7^ 0; cf. [20]. One may treat the general particular solution of the form (19. 2p with 
the coefficients (I9.6p - (l9.8p and (19. lip as an example of application of yet unknown "nonlinear" 
superposition principle for the Schrodinger equation under consideration for two particular solutions 
of a similar form with ci 7^ 0, C2 = and ci = 0, C2 7^ 0. 

It is worth noting that function (19.20 with the coefficients given by (I9.4p - (l9.1ip does also satisfy 
the following linear Schrodinger equation 

^^ = -a (t) ^ + & it) x'^P - I [c it) ^^ + dit)^)+-^^^, s> 0. (9.12) 
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Then a more general solution of this equation can be obtained by the superposition principle as 
follows 

/oo 
Kh{x,y,t) x{y) dy, (9.13) 
-oo 

where x is an arbitrary function such that the integral converges and one can interchange differen- 
tiation and integration. Solution of the Cauchy initial value problem simply requires an inversion 
of the integral 



V'(a;,0)= / Kf,{x,y,0) xiv) dy (9.14) 

J — oo 

that is 

X (y) = r Kl (x, 0) ^ (x, 0) dx, (9.15) 

say, by the inverse of the Fourier transform. Thus our equations (I9.13P and (19.151) solve the initial 
value problem for the above linear Schrodinger equation (19.121) as a double integral with the help 
of the kernel Kh (x, t) that is regular at t = 0, when /i (0) = ci 7^ 0. 

On the other hand, 

K„{x,y,t)= / Gh{x,z,t) Kh{z,y,0) dz, (9.16) 



where Gh {x, y, t) is the Green function, which can be obtain from our solution (19. 2p in the limit 
Ci ^ with a proper normalization as in the propagator (11.51) . Therefore, substitution of (19.151) 
into (I9.14P gives the traditional single integral form of the solution in terms of the Green function 

/•oo 

^(x,t) = / Gh{x,y,t) ^(y,0) dy (9.17) 



by ^M). 

The case of the new Hamiltonian, corresponding to (12.30 . is similar. The general solution of 
characteristic equation (I2.14p is given by 

H = C2fi2{t)+C3fi3{t) (9.18) 
= cos t (c2 sinh t — cosh t) + sin t (c2 cosh t + sinh t) 

= -Le*(c2sin(t + ^)-C3Cos(t + £ 

--^ e^* cos + £ + C3 sin (t + 

and /i' = 2 cosht (c2 cost + C3 sint) with /i (0) = — C3, fi' (0) = 2c2. The first three coefficients of the 
quadratic form in the solution (19.20 are 

, , cos t (c2 cosh t — C3 sinh t) + sin t (c2 sinh t + C3 cosh t) . . 

it] = ( y . xy ) 

2 (cos t (c2 sinh t — C3 cosh t) + sin t (c2 cosh t + C3 sinh t)) ' 

Pit) = ^-T r?^^°^ , r -r-T, (9-20) 

cos if: (C2 smh t — C3 cosh tj + sm t (C2 cosh t + C3 smh t) 

^U) = 7(0) + C3/?^(0) (cos t sinh t + sin t cosh t) ^i) 

2 (cos t (c2 sinh t — C3 cosh t) + sin t (c2 cosh t + cs sinh t)) 
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and one can use formula (19. lip for the last coefficient. The corresponding elementary integral is 

dt sin t 



I 



+ C. 



{Acost + Bsintf A {Acost + B sint) 
The cases (I2.18P and (12.191) can be considered in a similar fashion. The results are 

fl = Cs/ig (t) + C4Ai4 (t) 

= sin t (c3 sinh t + C4 cosh t) — cos t (03 cosh t + C4 sinh t) , 
sin t (cs cosh t + C4 sinh t) + cos t (03 sinh t + C4 cosh t) 



a{t) 
Pit) 



l{t) = 7(0) + 



2 (sin t (c3 sinh t + C4 cosh t) — cos t (03 cosh t + C4 sinh t)) 

-C3/3 (0) 

sin t (c3 sinh t + C4 cosh t) — cos t (03 cosh t + C4 sinh t) ' 
€3/?^ (0) (sin t cosh t — cos t sinh t) 



2 (sin t (c3 sinh t + C4 cosh t) — cos t (03 cosh i + C4 sinh t)) 



and 



a(t) 
/3(t) 



= COS t (ci cosh t — C4 sinh t) + sin t (ci sinh t + C4 cosh t) . 

sinh t (ci cos t + C4 sin t) + cosh t (ci sin t — C4 cos t) 
2 (sinht (ci sint — C4 cost) + cosht (ci cost + C4 sint)) 

ci/j(0) 

sinh if: (ci sin t — C4 cos t) + cosh t (ci cos t + C4 sin t) ' 

(0) (cos t sinh t — sin t cosh t) 



lit) = 7(0) + 



9.22) 
9.23) 

9.24) 
9.25) 
9.26) 

9.27) 

9.28) 
9.29) 
9.30) 



2 (sinh if (ci sin t — C4 cos t) + cosh t (ci cos t + C4 sin t) ) ' 

respectively. One can use once again formula (19. lip for the last coefficient. We leave further details 
to the reader. 

10. A Note on The Ill-Posedness of The Schrodinger Equations 

The same method shows that the joint solution of the both linear and nonlinear Schrodinger 
equations (19.120 and (19. ip . respectively, corresponding to the initial data 

i^\t=o = Se{x-y) = -jL= exp ~ j , s > 0, (10.1) 

has the form 

ij = Ge {x, y, t) = ^ ^i{aAt)x^+pAt)-y+^At)y'+-s(t)) (iq_2) 

V^/^e (t) 

with the characteristic function fi^ (t) = 27t [sfii (t) + /ig (t)) . The coefficients of the quadratic form 
are given by 

cos t (e sinh t + cosh t) — sin t (e cosh t + sinh t) 

(t) = 7 7 7 TT-, (10.3) 

2 (cos t [e cosh t + smh t) + sm t [e smh t + cosh t) ) 
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1 

COS t {e cosh t + sinh t) + sin t {e sinh t + cosh t) ' 
cos t cosh t + sin t sinh t 
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(10.4) 
(10.5) 



2 (cos t (e cosh t + sinh t) + sin t (e sinh t + cosh t) ) 

We simply choose ci = 27r£: > 0, C2 = 27r and e**^ = 1/v^ and the initial data a (0) = 7(0) = 
—P (0) /2 = 1/ (2e) in a general solution fl9.6p - fl9.8l) . The case e = 0, t > corresponds to the 
original propagator (11. 5p . while e > 0, t = gives the delta sequence (110. ip . 

If h = hs{t) = ( A/27r) /i^ = 2 A cos t {e sinh t + cosh t) , then 



X (e/ii(t)+/i2(t)) 



l~s 



A-s 



1 - s 



when < s < 1, 
when s = 1 



(10.6) 



with Kg (0) = provided e > 0. 

In this example, the initial data "^/'l^^o ~ Ge{x,y,0) = 5^{x — y) converge to the Dirac delta 
function 5 [x — y) as e — >■ 0"*" in the distributional sense [17], [58], [M], [66] 



On the other hand, 



lim 



lim 



Ge {x, y, 0) {y) dy = ip (x) . 



Ge{x,y,t) ip{y) dy 



[10.7) 



[10.8) 



Go{x,y,t) ip{y) dy 



with t > 0. When s = 1 the solution ip = Ge {x,y,t) , t > does not have a limit because of 
divergence of the logarithmic phase factor it) as e 0^. See also Refs. [3] and [33] on the 
ill-posedness of some canonical dispersive equations. 

The second case, corresponding to (12.31) . is similar. One can choose /i^ (t) = 2ti (/ig if) — efi^ (t)) 
and obtain 

cos t (cosh t + e sinh t) + sin t (sinh t — e cosh t) 



as (t) 

it) 



2 (cos t (sinh t + e cosh t) + sin t (cosh t — e sinh t) ) 

1 

cos t (sinh t + e cosh t) + sin t (cosh t — e sinh t) ' 
cos t cosh t — sin t sinh t 



2 (cost (sinht + ecosht) + sint (cosht — esmht)) 
If he (t) = (A/27r) /i^ = 2A cosht (cost — e sint) , then 

A (/i2 (t) - £/i3 (^)) 



(10.9) 
(10.10) 
(10.11) 



He (t) 



\1 — S 



(2vr)^ 



'2^ 



In 



when < s < 1, 
when s = 1 



(10.12) 
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and (0) = when e > 0. Formulas fll0.2p and fll0.9p - fll0.12l) describe actual (nonlinear) evolution 
for initial data as in fllO.ip . One can observe once again a discontinuity with respect to these initial 
data as e ^ 0^. 



The cases (12. 180 and (12.190 are as follows. One gets fi^ (t) = 27r (/i4 (t) — efi^ (t)) 

sin t (sinh t — e cosh t) + cos t (cosh t — e sinh t) 



and hs (t) 



as (t) 

(A/27r)/i; = 
^^e {t) = 



2 (sin t (cosh t — e sinh t) — cos t (sinh t — e cosh t) ) 

1 

sin t (cosh t — e sinh t) — cos t (sinh t — e cosh t) ' 
cos t cosh t — sin t sinh t 
2 (sin if: (cosh t — e sinh t) — cos t (sinh t — e cosh t) ) 

2A sin t (sinh t — e cosh t) , 



A (/i4 (t) - e/ig (t)) 



(2^)' 

27r 



^4 W 



1 - S 



when < s < 1, 



when s = 1 



with (0) = 0, e > in the case fIXTSD . Also /x^ (t) = 27r (/i4 (t) + eyU^ (t)) , 

sinh t (sin t + e cos t) — cosh t (cos t — e sin t) 



and (t) 



tte (t) 
it) 

(A/27r) = 

f^e it) = 



2 (sinh t (cos t — e sin t) — cosh t (sin t + e cos t) ) ' 

1 

sinh t (cos t — e sin t) — cosh t (sin t + e cos t) ' 
cos t cosh t + sin t sinh t 
2 (sinh t (cos t — e sin t) — cosh t (sin t + e cos t) ) 

2A sinh t (sin t + e cos t) , 



(2 



27r 



/i4 W 



1 - S 



^1 (i^) 



when < s < 1, 



when s = 1 



;i0.13) 
;i0.14) 

;io.i5) 



^10.16) 



(10.17) 
(10.18) 
(10.19) 



[10.20) 



with Kg (0) = 0, e > in the case (12.191) . We leave the details to the reader. 
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11. Appendix A. Fundamental Solutions of The Characteristic Equations 
We denote 

Ui = cost, M2 = sint, fi=cosht, f2 = sinht (H-l) 
such that u[ = —U2, u'2 = Ui, v[ = V2, v'2 = fi and study differential equations satisfied by the 
following set of the Wronskians of trigonometric and hyperbolic functions 

{W K, t;/3)},,^=i,2 = ^1) ' ^ ^2) , W {U2, Vi) , W {U2, V2)} . (11.2) 

Let us take, for example, 

y = W {ui,vi) = U1V2 + U2V1. (11-3) 

Then 

y' = 2uivi, y" = 2uiV2 - 2u2Vi (11-4) 

and 

y" - Ty' + Aay = (4cr + 2) uiV2 + {Aa - 2) U2V1 - 2tuiVi = 0. (11.5) 
The last equation is satisfied when a = 1/2, r = 2v2/vi and cr = —1/2, r = —2u2/ui. All other 
cases are similar and the results are presented in Table 1. 

Our calculations reveal the following identities 

W" {Ui.Vi) = -2W {U2,V2) , W" {Ui,V2) = -2W {U2,V{) , (11.6) 

W" (m2, V,) = 2W (Ml, V2) , W" {U2, V2) = 2W (m, V,) , 

for the Wronskians under consideration. This implies that the set of Wronskians flll.2p provides 
the fundamental solutions of the fourth order differential equation 

iy(4)+4iy = (11.7) 

with constant coefficients. The corresponding characteristic equation, + 4 = 0, has four roots, 
Ai = 1 + i, A2 = 1 — ^, A3 = — 1 + z, A4 = — 1 — z, and the fundamental solution set is given by 

W^/3}a,/3=l,2 = {"1^1' "1^2, U2V1, U2V2} ■ (11.8) 

These solutions of the bi-harmonic equation flll.Tp are even or odd functions of time. They do 
not satisfy our second order characteristic equations. For example, let Wi = U1V2 = costsinht and 
Wi = U2V1 = sin t cosh t. Then, by a direct calculation, 

sinh t 

L (wi) = w'l + 2tant w[ - 2wi = -2 , (11.9) 

cost 

,, , sinht 

L {W2) = Wo + 2tant — 2^2 = 2 . 

cost 

Thus, separately, these solutions of (111.71) satisfy nonhomogeneous characteristic equations. But 
together, 

^inh / sinh t 

L (y,) =L{w^ + W2) = L (w^) + L (102) = -2—— + 2 = 0. (11.10) 

cos t COS t 

A similar property holds for all other solutions of the characteristic equations from Table 1. 



28 



RICARDO CORDERO-SOTO AND SERGEI K. SUSLOV 

Table 1. Fundamental solutions of the characteristic equations. 



(^naidcteiibtic equation y — ly -\- ^oy — u 


runQdnientdi solution set ykSi^k 


u" + U = Q 

{a — i/4, r — uj 


ui = cos t, U2 = sin t 

\Ui — —U2, U2 — U\) 


v" -v = Q 

(rr — —1 /A T — r\\ 


Vl — cosh = sinh t 

— V2i U2 — i>i) 


l/" + 2tant|/'-2y = 
(a = -1/2, r = -2m2/mi) 


y^ = W {Ui,Vi) = U1V2 + U2V1 

y2 = W {ui, V2) = uiVi + U2V2 


?/"-2coti y'-2y^0 

{a = -1/2, r = 2m2/mi) 


ys^W {u2, V2) = 1*2^1 - U1V2 
y4 = W{u2,vi)= U2V2 - uivi 


?/"-2tanht y' + 2?/ = 
[a = 1/2, r = 2t;2/t;i) 


yi = W {ui, Vl) = U1V2 + U2V1 

y^ = W {u2,Vi) = U2V2 - UiVi 


- 2cotht y' + 2y = 
{a = 1/2. - = Orjr,) 


y2 = W {Ui,V2) = + U2V2 
y?, = ("2- ''2) = »2''l - »l''2 



In order to obtain the fundamental solutions in an algebraic manner, we denote 

L2 



Li = — + 2^--2, 
dt"^ Ui dt ' 

L3 = ^-2^- + 2, 
dt"^ Vl dt ' 



d^ ^ui d 

2— 2 

dt"^ U2 dt ' 

d^ Vl d 
2— h 2 

dt"^ V2 dt 



:ii.ii) 



and compute the actions of these second order linear differential operators L^. on the four basis 
vectors {uaVfj}^ ^ , namely, {uaVp) . The results are presented in Table 2. 

Table 2. Construction of the fundamental solutions. 



Linear operators 


UiVi 


U1V2 


U2V1 


U2V2 


d 

It 


UIV2 - U2VI 


UiVi - U2V2 


UiVi + U2V2 


U1V2 + U2V1 


d^ 
dt^ 


-2U2V2 


—2u2Vi 


2uiV2 


2uiVi 


Li 


-2^ 

Ui 


V2 
-2 — 

Ui 


2^ 

Ui 


2^ 

Ul 


L2 


V2 
-2 — 

U2 


v^ 

-2 — 

U2 


-2^ 

U2 


V2 
-2 — 

U2 


Lz 


2^ 

Vl 


-2^ 

Vl 


2^ 

Vl 


2^ 

Vl 


L4 


2^ 

V2 


-2^ 

V2 


-2^ 

V2 


-2^ 

V2 



Therefore 



Li {uiVi + U2V2) = Li {uiV2 + U2V1) 

= L2 {U2V2 — UiVi) = L2 {U2VI — U1V2) 

= L3 (iiii;2 + 1^21^1) = L3 (U2V2 - UiVi) 

= L4 (lii'Ui + W2'y2) = L4 {U2V1 — U1V2) = 



;il.l2) 
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as has been stated in Table 1. 

All our characteristic equations in this paper obey certain periodicity properties. For instance, 
equations 

y" + 2ia.nt y' -2y = Q (11.13) 

and 

y" -2coit y' -2y = Q (11-14) 

are invariant under the shifts t ^ t ±tt and interchange one into another when t — > t ± 7r/2. Since 
only two solutions of a linear second order differential equation may be linearly independent, the 
corresponding fundamental solutions satisfy the following relations 

yi (t ± tt) \ _ / cosliTT isinhvr \ / yi{t) 
y2it±n) J \^ isinhvr coshvr J\y2it) 

and 

f y,{t± 7r/2) \__f sinh {rr/2) ± cosh {n/2) \ f y^it) \ 

[ y2{t±n/2) J~ [ ±cosh(7r/2) sinh (7r/2) j[y,it) J ' ^^'"'^^ 

respectively. Two other characteristic equations have pure imaginary periods. We leave the details 
to the reader. 



;ii.i5) 



12. Appendix B. On A Transformation of The Quantum Hamiltonians 

Our definition of the creation and annihilation operators given by (11. 3p implies the following 
operator identities 

a;2 = i [aa^ + a^a) - ^ + (a^Y) , (12-1) 
^ = (aat + ata) - \ {a' + (af) , (12.2) 

2x^ + 1 = -a' + {a^Y (12.3) 

(and vise versa), which allows us to transform the time-dependent Schrodinger equation (12. ip into a 
Hamiltonian form (11.11) . where the Hamiltonian is written in terms of the creation and annihilation 
operators as follows 

H = ^{a{t) + b{t)) {aa^ + a^a) (12.4) 

+ i (a (t) -b{t) + 2id{t)) + ^ (a (t) -h{t)- 2id{t)) {a^f , 

when c = 2d. This helps to transform the Hamiltonians of modified oscillators under consideration 
into different equivalent forms, which are used in the paper. 

The trigonometric cases (12.20 and (12.180 results in the Hamiltonians (11.21) and (17.181) with n = 1, 
respectively. In the first hyperbolic case (12.31) one gets 

1 1 1 2 

H = - cosh {2t) [aa) + a) a) + -{I - i sinh (2t)) + -{I + i sinh (2t)) [a)) , (12.5) 

where 

1 ± i sinh {2t) = cosh (2t) e±^'^'-^*'^'^(2r)^ ^ _ 1 ^^^^ ^2t) , (12.6) 

2 
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which imphes the Schrodinger equation (ll.6p - fll.7p . The second hyperbohc case (12. 191) is similar. 
We leave the details to the reader. 



13. Appendix C. On A Hamiltonian Structure of The Characteristic Equations 

The Hamilton equations of classical mechanics [35] , 

. dH . dH 

g = ^, P=-^, (13.1) 
op oq 

with a general quadratic Hamiltonian 

H = a{t)p^ + h{t)q^ + 2d{t)pq (13.2) 

are 

q = 2ap + 2dq, p = -2bq - 2dp. (13.3) 

We denote, as is customary, differentiation with respect to time by placing a dot above the canonical 
variables p and q. Elimination of the generalized momentum p from this system results in the second 
order equation with respect to the generalized coordinate 

g--g + 4la6-rf2 + ^l--3llg = 0. (13.4) 




It coincides with the characteristic equation fl2.8l) - fl2.9l) with c = 2d. Our choice of the coefficients 
fl2.2p - fl2.3l) and fl2. 180 - 02. 190 in the classical Hamiltonian (113.20 corresponds to the following models 
of modified classical oscillators 

g + 2tantg-2g = 0, (13.5) 
g-2tanhtg + 2g = 0, (13.6) 

q-2cott q-2q = 0, (13.7) 
g-2cothtg + 2g = 0, (13.8) 

respectively; see Appendix A for their fundamental solutions. 

The standard quantization of the classical integrable systems under consideration, namely, 

d 

q ^ X, p-^i^^T^, \x,p] = xp — px = i (13.9) 
ox 



and 



dip 

H ap^ + hx^ + d {px + xp) , z— = Hil), (13.10) 
leads to the quantum exactly solvable models of modified oscillators discussed in this paper. 

Another example is a damped oscillator with a = (a;o/2) e~^'^*, b = (a;o/2)e^'^* and c = d = 0. 
The classical equation 

q + 2Xq + ujlq = (13.11) 
describes damped oscillations [35j. The corresponding quantum propagator has the form (12.40 with 

^ = ^e'^*smut, uj^ = iul - X'^ > (13.12) 
u 
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and 



, , LJ COS out — X sin out 
a(t) = e^^*, 13.13 



The Schrodinger equation 



^t) = ^e'\ (13.14) 

ujq sm uji 

u"^ — Un sin^ ut 

7 t = 5 ^ -. 13.15 

zujq sm ujt [LJ cos LJt — A sm ut) 



At 



describes the hnear oscillator with a variable unit of length x —>■ xe . See [21] for more details 
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